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I. INTRODUCTION 


Quantum topological field theories, such as the Chern-Simons and BF models, are gauge- 
invariant and metric-independent theories. Due to the latter characteristic, they can be used 


for the study of knot and link invariants in three dimensional {3d) spaces 0,0,1 , and for 
the study of generalizations of these invariants appropriated to the dimension of the base 
manifold. 

This relationship between topological theories and knot invariants subsists even at the 
classical level, as has been shown in the Abelian case and also in the non-Abelian case 
0 by means of a perturbative study of the classical equations of motion of topological 
theories coupled to external point-particles carrying non-Abelian charge (Wong particles 
[^). The method presented in these references rests upon the fact that the classical action 
of the theory must retain its diffeomorphism-invariant character when it is evaluated on- 
shell. Hence, the on-shell action of the Chern-Simons (or BF) theories coupled in a suitable 
manner to particles should yield link invariants in 3d, just as the vacuum expectation value 
of the Wilson-Loop does within the quantum held approach. This result concerning the 
classical perturbative treatment can be rigorously proven and generalized to situations A^ere 
the symmetry group is other than the group of diffeomorphisms of the base manifold 


Q. 


The purpose o 
in references 


this article, within the context of the ideas and procedures developed 
^,0], is to apply the classical perturbative method outlined above to 2d 
Euclidean Yang-Mills theory, which yields and example of how the method works when the 
symmetry group is smaller than the full group of diffeomorphism. As it is well known, 2d 


Yang-Mills is invariant under diffeomorphisms that preserve areas jalSlllOl, which reflects 
in the fact that the Wilson loop exhibits an exponential dependence on the areas of the loop 
ll|. 2d Yang-Mills theories have been studied in many contexts. For instance, they describe 
closed topological strings with bounded states, and q-deformed 2d Yan g M ills theories in 


Riemann surfaces give topological invariants in one dimension higher 


QQ. 


Furthermore, 


under certain conditions. Euclidean 2d Yang-Mills theory mimics the conhning phase of 
QCD to a good degree of accuracy QQ. 

The p^er is organized as follows. In Section II we apply the general method developed 
in 4 S Is! to 2d Yang-Mills theory coupled to external Wong-particles. After presenting 
the model we discuss a perturbative scheme for solving the classical equations of motion. 
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and calculate the on-sliell action to the hrst two orders in the coupling constant. In section 
EH we hnd a geometric interpretation of the invariants. We also discuss the consistency 
conditions that gauge invariance imposes over the classical equations of motion, and their 
relationship with the dehnition of the surface-invariants. To this end, we perform a gauge 
transformation that amounts to attaching a bundle of straight lines to each point of space. 
In this gauge, area-invariants become neatly expressed in terms of the areas formed by the 
cross product among the tangent vectors to these “hbers” and those tangent to the curves. 
Some hnal comments close this section. 


II. YANG-MILLS THEORY IN 2 DIMENSIONS COUPLED TO WONG PARTI¬ 
CLES 


A. The model and the general method 


Our starting point will be the action 


S — SyM + Sint — 


^ [ d^xTr{F^nFn + ^2 f drTr{K,g-\T)D, 

^9 J J-Y 


( 1 ) 


where the last term corresponds to the interaction of n non-dynamical Wong particles (i.e.. 


classical particles with cromo-electrical charge) 0. O' is the gauge coupling constant of the 
Yang-Mills held. Since we are dealing with the Euclidean theory, there is no distinction 
between “Greek” (/r, z/, etc.) sub or super-scripts, so we shall use both indistinctly. We will 
use the following conventions for the — 1 generators T® of the su{N) algebra and for the 
gauge held 


Tr{T^T^) = 

( 2 ) 

j'^a _ j^abcrjnc 

( 3 ) 


( 4 ) 

Ai = A^(Zi(T))z^(T). 

( 5 ) 


As dynamical variables we take the potentials A“ and the SU{N) matrices in the fun¬ 
damental representation gtir), associated with the internal degrees of freedom of the Wong 
particles |^. Instead, the trajectories 2 ;f(r) of the particles are given. The curves drawn 
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by the particles are then taken as closed curves in the Euclidean plane, surrounding surfaces 
Sj whose “area-invariants” we are going to study. 

In equation ([Q) we use the covariant derivative along the world-line of the i-particle: 
DrSiir) = gi + Aj(r)( 7 i(r). Also, it appears Ki = which is a constant element of the 

algebra related to the initial value of the cromo-electric charge liir) through 

hir) = gi{T)Kig^\T). ( 6 ) 


Finally, the Yang-Mills held tensor is dehned as usual by 


FfMu — -f [A^, Ay], 


( 7 ) 


The hrst term in the action (i.e., the Yang-Mills one) is invariant under area-preservi^ 
diffeomorphisms. The second one is invariant against arbitrary diffeomorphisms 3, l2|- 
Hence, the whole action is invariant under area-preserving diffeomorphisms. Also, it is 
invariant under gauge transformations 


A„ A'J = 


Ki ^ K ‘ = 


9i 

h 


9i = 


j\ I _ 


K^, 

gi, 


( 8 ) 

( 9 ) 

( 10 ) 

( 11 ) 


In can be seen that the interaction term of the action m may be put in the form 
Sint = drzi^r) ^r[Kig~^df,gi\ + j (fx6^‘^\x - Zi{T))Tr[Ii{T)Af, (zi^r))'^ , ( 12 ) 

hence, varying the action (CD with respect to the dynamical variable A“ we obtain 


= Ar = / dTz''{T)hiT)6^‘^\x-z,{T)), 

i=l dzi 


(13) 


where we have set A = g'^. 

To write down the equations of motion for the internal variables giix) we follow the 
procedure given in Take a parametrization of the group elements 


9i = 9i{^i) = 


(14) 
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and perform variations of Sint with respect to the — 1 independent parameters 
This leads ns to the Enler-Lagrange eqnations 


dL 


d / dL 




dCfir) dT\dif{T) 


= 0 , 


where we have dehned 


i = '^Tr{Kig^ ‘(T)£>^9j(T)). 


( 15 ) 


( 16 ) 


It can be seen that eqnations m are eqnivalent to the gange-covariant conservation of 


the non-Abelian charge of each particle along its world line 
by taking the covariant derivative on both sides of eqnation m: 


Ifij l (this conservation-law arises 


DtU — a + [Aj, /j] — 0, 


(17) 


whose solntion is 

UT) = U,{T)mUr\T). ( 18 ) 

Here t/j(r) is the time ordered exponential of the gange held along the curve 7 * 

Uiir) = Texp (- f Aiir') dr' ). (19) 

Jo 

From and (HHD we obtain gtir) = f/j(r) 5 'j( 0 ), which in turn implies Dngiir) = 0. There¬ 
fore, we hnd that the interaction term Sint of the action vanishes when it is evaluated 
on-shell. Hence, we only have to consider the Yang-Mills action evaluated on the equations 
of motion. 

Plugging m. and (USD into m yields an equation for the gauge potentials in terms 
of the curves 7 * j^. Inserting the solution to this equation into the action (HD one would 
hnally obtain a functional S'([ 7 ];A) that only depends on the curves 7 and the coupling 
constant A P| . Following the general arguments discussed in the introduction, this on-shell 
action should retain the invariance under area-preserving diffeomorphisms. Henceforth, the 
hnal expression for S'([ 7 ];A), calculated as explained above, constitutes an invariant under 
area-preserving diffeomorphisms associated to the curves 7 *. 


B. Perturbative expansion to the first two orders 

Since the equations to be solved are non-linear, we will use a perturbative method to deal 
with them, such as in reference a. A5 la., , wUI ..a, ,.at action on-aKell tnay be 
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written as a power series in the parameter A 


A 


Son—shell (|7.].a) = -5; A” Sl'’>|7j 


( 20 ) 


p=0 


In order to carry out the perturbative method, we dehne quantities: 

m = A"^A., 




R7 ^ 


( 21 ) 


The equations (HZD for the parameters, and their formal solution (uni), adopt the form 

dl-ir) 


dr 


+ ARr(T)r:(T) = 0 


( 22 ) 


and 


/“(t) = [Texp ( - a / i{i(r')dT')]“/nO). 


(23) 


respectively. We need to develop the ordered exponential in the right hand side of m- 

n „ 

i=l -Izi 

-A^j2 / ^“‘(7i)'5‘"’(^-‘-)A“‘(0) 

i=l •’Zi do 

^ e pZ pZ\ 

i=\ 1i ^0 ^0 


+(-A)'+‘^ I dz- fdz^,'... r"dzf/RZ'(zO 

I—I d-fi Jo do 

■ ■■RZ-'^(z,)S‘'^>(x-z)It’(0) 

; (24) 


Also, we expand the “new” helds in powers of A, 

OO 

“; = E(25) 

p=0 

and introduce them into ()24|1 to obtain the right hand of m order by order. Additionally, 
we have to work out the left hand side of m 

= DA" - + a^[A^ a"] + [ a „ - a^A^] + [a^, [a^ a"]] , ( 26 ) 
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where □ = (since we are in Euclidean space this is the same as the Laplacian V^). 

Choosing the Lorentz gauge = 0, and changing to the variables in (EHD we can write 
the equation of motion (USD as 

An< + A^2a^,dXfabc - + A'a^fea^^/cde/abc = AJ^, (27) 

where the right hand side of this equation is given by (1^ . with the substitution (ES) , as 
discussed. 

We are ready to evaluate the Yang-Mills action on-shell up to the hrst orders in A. It will 
be useful to dehne the ’’Abelian part” of the held tensor which allows to 

write the Yang-Mills action as 

-SvM = ^ / d^xTr{F^,Fn = ^J d^xTr{f^, + A[a^, a,]f 

= d‘^xTr{ff,^f^'' + 2Afi„y[af„a^] + A‘^[ai,,a,y][a>^,a’']). (28) 

From expression it is immediate to obtain the 0 — th order contribution of the on-shell 
action CD 

4m = ^ / iPxTr(FlSF'^<'») = ^j XxTr(U,D = j ^xTr(a,a>«> - 

(29) 

where is the solution to the 0 — fh-order equation of motion that results from (ESD 

n « 

□^M(O) _ JM(0) _ ^ i dz^h(^)(x - z) /“(O). (30) 

j=i 

Before showing the hnal expression for the 0 — fh order invariant, we hnd it convenient to 
work out the hrst order invariant to the same level of detail. From (EHp we have 

4i = j rf^lTr (/«/<“■<»)+/(»)K<»),a-<»)l) 

= y'<i=irr(-2a<“>na<‘' + /l»)[a'‘<»),a-<»)|), (31) 

where we have thrown away boundary terms and used the Lorentz gauge. From (EHD , we 
have that the hrst order potential obeys the equation 

( 32 ) 
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where can be easily obtained from (El- 

At this point it is worth noticing that the general strnctnre of the p — fh-order eqnation 
of motion is, indeed, the same as that of the hrst two ones: the Laplacian of the p — fh-order 
potential is given in terms of previons orders potentials, which were already solved as 
fnnctions of the cnrves 7 j. Therefore, the pertnrbative method can be recnrsively applied, 
and the on-shell action can be obtained order by order. 

Tnrning back to the hrst contribntions, we see that it remains to solve the 0 — fh and 
hrst orders eqnations (pn|l and (jS 21 ), and to snbstitnte these resnlts into and (jS!}. To 
this end, we hnd it convenient to introdnce loop-coordinates |l^ 

Jo 

rzi rZn -1 

X / dz!^^... - z)5^^\x 2 - zi)5^^\xs - Z 2 )... M^\xn - z^-i), 

Jo Jo 

(33) 


that obey the diherential constraints 

^ ^ (-6{Xi - Xi_i) + 6{Xi - Xj+i)) ^ 3 ^^ 




where both xq and Xn+i are taken as the origin of the loop. Loop-coordinates also obey the 
algebraic constraints 


+ ' 'rj^Pk(Pl--Xn) — JlMl'''Mfc + l'''Mn 

Pk 


(35) 


The snm in the right hand side of the last eqnation is made over all the permntations of 
the indices p that preserve the order on the snbset pi - ■ ■ pk and on the remaining snbset 


hfc+l ' ' ' dn- 

We shall nse a generalized Einstein convention 




Y, I A^^B^-'^y-d^x = Y I A^ix)By^-ix,y...)d^ 


X. 


(36) 


II II 

A bar over a“continnons index” breaks the Einstein convention 


A^.^yB^^'^y = Y [ Ay.^yBy^’^yj^x. ( 37 ) 

For the sake of simplicity, we restrict onrselves to consider SU{2) as gange gronp. We 
shall nse arrows to denote iso-vectors, and employ the ’’dot” and ’’cross” prodnct notation. 
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Using these tools, and can be written down as 

sfl, = 2 j SxdJ°^ ■ 

= - 24 »).na<,»l, 


and 


( 38 ) 


(39) 


s?M = 4 "i • + - (42 X 42 ). 

while the 0 — th and 1 — st order currents are given by 

J{0)iix _ E T^h 

i 

j(lUx _ ^Tl^xuy^ifi) ^ 
i 

In turn, equations (EDI) and (El for the zero and hrst orders can be written as 

i 

and 

(b = (0) X - 2#^ x ^ x k. (42) 

i 

The solution of dH is given by 


^(0)/ix 


= G.., = — 5] 4" In |i - !/| t, 


(43) 


In this equation, is the Green function of the Laplacian: ‘^6^‘^\x — y) = ^ In \x — y\ = 


Gx,y 


Substituting El in El we obtain the area-invariant corresponding to the 0 — th order 


S- 


( 0 ) 

YM 


On—Shell 


- Elt ■ tl - 31 h") = 2 Elt ■ tl 3j). (44) 

ij i,j 

where we have dehned the quantity J{pfi,^j). It is worth noticing that this expression 
also corresponds to the result that we would obtain if we had considered the Maxwell theory 
coupled with Abelian particles. In section mTl we shall interpret the meaning of this invariant. 

Introducing (El and El into El we can also obtain the hrst order Yang-Mills on-shell 
action in terms of the curves The result is 




( 1 ) 

YM 


On-Shell V 2 n 


Yllhry-!,) ■ h] b 

i,j,k ^ 


_rj:jlXlrpUX2rpUy | 111 \x — 1/| 

^nxtxyj.yx^j.ux 2 In |y - Xi | 111 |x - Xa | } • (45) 
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III. INTERPRETATION OF THE INVARIANTS. CONSISTENCY CONDITIONS 


A. Interpretation of the 0 — th order invariant 


Since the isovectors are arbitrary we conclude, from (m, that the quantities 

jhi, ij) = 111 1 ^ - 2/1 = -^ / 2 / 2 ~ 

must be area invariants. In order to interpret them we introduce the 0-form 

F,,s = F(y,S)^ [ £x6^^\x-y), (47) 

JT. 

with support on the 2-dimensional surface S bounded by 7 . It can be readily seen that 

e^^d,yF{y,T.)= [ dx^d^^\x-y). (48) 

J-l=dT, 

Substituting the loop coordinates in dH with the help of the previous expression we 
obtain 




- 7 ^/ d'^x [ d‘^yeyad°‘yFy^j:An\x-y\e^^dpsiFs,j:j 

-TT I I In 1^ - y\Fy,^i Fx,^j 



(49) 


Expression dH measures the common area of the surfaces Sj and Sj bounded by the 
closed curves 7 * and 7 ^ respectively. Clearly, this quantity is invariant under area-preserving 
diffeomorphisms, as expected. It is interesting to notice that Jhi,lj) formally analogous 
to the Gauss linking number of two curves in [^, which is just the 0 — f/i order link 
invariant obtained when the perturbative procedure discussed here is applied to the classical 
Chern-Simons theory 

Unfortunately, the “trick” of substituting the loop coordinate using ()48|1 and in¬ 
tegrating by parts does not suffice to yield a geometric interpretation of the hrst order 
area-invariant (gSl), since it is not possible to get rid of the logarithm functions that appear 
in the expression for that invariant. Later, in subsection 1111 Cl we shall present a different 
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method that provides a geometric interpretation of the hrst order invariant, and hopefully 
of the higher order ones too. Before, we study the consistence of the 0 — th and hrst orders 
equations of motion, since that method relies on the existence of their solutions. 

B. Consistency conditions 

The consistency of the 0 — fh order equation of motion dn} is immediate: taking the 
divergence of both sides of the equation, and using our gauge condition, we obtain 

= Y, d.rrh = 0, (50) 

i 

since we are taking the trajectories of the Wong particles as closed curves. 

Regarding the hrst order equation dSl , we rewrite it as 

^ _ 2-(o)mx ^ ^ 

Taking the divergence on both sides of m and using the gauge condition we obtain 

X x x = 0. (52) 

Using expressions iOD, (gH) and dn for jdU®, □^( 0 )^ 3 : respectively, and taking 

into account the diherential constraint dMD we have, after some calculations 

(/;• X - x,{0)) rr In |x - y\ = 0, (53) 

ij 

where Xj(0) is the starting point of the curve 7 *. If the curves do not intersect each other, 
and their iso-vectors are independent, equation d^ leads to the condition 

•^(T,7i) = 0, (54) 

for every i,j. Hence, the hrst order equation of motion is consistent, and consequently the 
hrst order contribution to the action is meaningful (for arbitrary values of the iso-vectors 
li) just when the 0 — th order invariant J( 7 j, gj) vanishes for every i,j. 

It is interesting to notice that this result resembles what occurs in the study of link 
invariants in through this method P|: the next order invariant is meaningful whenever 
the preceding one vanishes. In turn, this coincides with what happens with the Milnor 
link invariants [l^, which suggests that the chain of area invariants that one would obtain 
following this method could be seen as a kind of “projection” of the Milnor invariants, from 
links in and diffeomorphisms, to surfaces in R^ and area preserving diffeomorphisms. 
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C. Axial gauge and bundle of parallel curves 


Let us assign a straight line to every point x of space. For definiteness, take these lines 
starting at the spatial inhnity, running parallel to the vertical direction from bottom to top 
until each one reaches its associated point x. The line associated to x is denoted by The 
one-index loop-coordinate T^{x, A^) of A^ is 

T>^{x,Xy)= [ dz>^6^^\x-z). (55) 

J\y 

From the very structure of the perturbative equations of motion (see comment after 
equation (1221)) Lorentz gauge = 0, it is clear that the Green 

function Gx,y enters in the solution of these equations only through its gradient dxyGx,y 
Now, it is easy to see that the substitution 

d^^G,,y T>^{x,xy), (56) 

induces an Abelian gauge transformation on the potential 

a>^ip) ( 57 ) 

where "0 is a certain function depending on the bundle of curves A. Moreover, it can be 
shown that if 77 ^ is a vector in the direction of A, = 0. Hence, the replacement (15611 

amounts to changing from Lorentz to axial gauge [actually, since space is Euclidean, Lorentz 
and Coulomb gauges are identical, as well as axial and temporal gauges]. Finally, it can 
be shown that (at least) the 0 — fh and hrst order contributions to the action on shell are 
invariant under Abelian gauge transformations, provided that their respective equations of 
motion are consistent. Thus we arrive to the following result. Although the calculations 
are done in the Lorentz gauge, we are allowed to perform an Abelian gauge transformation 
whenever the consistence conditions holds (it seems that this result holds to every order, 
but we do not have a proof yet). Changing to the axial gauge will allow us to interpret the 
Q — th and hrst order invariants. Let us begin with the former. 

Starting from expression (HED one has 

= -^e^^do.xF^^xT^y\xi\x-y\ 

= e^^F^^xT^yd^xGx,y 

= e^“Fs,.7X"^„(x,A"). (58) 
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Now, observe that the integration in x must be restricted to the area surrounded by 
curve i, due to the factor in the integrand. Furthermore, for each x, the presence of the 
factor Tq(x, X^) ensures that there will be contributions only when this point x is connected 
through a vertical line with a point y lying in the other curve j, in such a manner that the 
inhnitesimal displacements along the vertical path (dz") and along the closed curve j {dy^) 
form a non-degenerated parallelogram (see FIG.P). Under these conditions, the contribution 
associated to the point x is just the area of that inhnitesimal parallelogram: e^°‘dy^dz'^. 



FIG. 1: Configuration with a non-vanishing 0 — th order contribution 

On the other hand, for each x belonging to Sj, but not to the intersection (see x’ in FIG. 

there would be two contributions (more generally, an even number of them, which could 
be none as in the case of x”, shown in FIG. Q) because the parallel curves intersect twice 
the closed trajectory j. But in these cases, the net contribution of the “upper side” of curve 
j (i.e. the piece from A to B in clockwise sense) plus the “down side” (the part from B to 
A in counterclockwise sense) cancel each other. 

Hence, we have to sum up as many inhnitesimal parallelograms as necessary to cover the 
area of intersection, and we recover our previous interpretation of the 0 — th order invariant. 
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D. Interpretation of the 1 — st order invariant and final remarks. 

Our starting point will be expression (03) for the Yang-Mills action on-sliell up to first 
order, expressed entirely in terms of the closed curves. Since the coefficient related to the 
iso-vectors /j is completely antisymmetric in i, j, k, and these iso-vectors are arbitrary, we 
conclude that the following quantity is invariant: 

Qrpr/‘‘%‘'’lnlx-x,llnlx-X2ld„lnlx-pl 

In \y - III In li - I 2 I) . (59) 

Here, square brackets denote anti-symmetrization. 

Consider the first term in (El. Making the substitution ()48|1 . and after some calculations 
it results 


In |x - xil In |x - X 2 \d^^ In |x - y\ 

= d^FY.^y In |x - Xi| In \x - In |x - y\ 

= In |x - Xi\d0x2 In \x - X2\dpydyx In |x - y\ 

= -{27rfe^'"Fj:^xiFT.jX2Fj:^ydax^Gx,x^dyx2Gx,x2di3ydyxGx,y. (60) 

Making the substitution (El, i.e., changing to axial gauge, the last line of (El can be 
written down as 


- X^^)Ty{x, \^^)dpyTy{x, A"). (61) 

This quantity vanishes, because the straight lines are parallel, which implies that the 
cross product of their tangent vectors is zero. 

It remains to compute the second term in (El. Changing to the axial gauge by means 
of substitution El, we obtain 

y(7„7i,7») = + [tjk], (62) 

where [ijk] means that the whole expression must be anti-symmetrized in the indices inside 
the bracket. This symmetry property of El imply that the minimal number of curves that 
will produce a non-trivial result is three. 
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To interpret this expression, let us then consider a configuration of three curves as in FIG. 
121 This configuration was chosen taking into account that, as we already know, 1 ^( 7 *, 7 j, 7 ^) 
is meaningful whenever (gni) vanishes for every pair of curves. In this configuration, the 
absolute values of the areas of intersection Aj, i = 1, 2, 3,4 are the same, but their signs are 
such that the above condition is fulfilled. It should be noticed that, given the configuration 
of FIG. 121 only the first term in contributes. The points xi and X 2 must belong to 



FIG. 2: Example of a “linking-area” configuration detected by F( 7 j, 7 j, 7 fc) 

the surfaces and respectively, as indicated by the surface-functions and Fy,^x 2 - 

Also, we have now two bundles of open straight parallel lines: the “fibers” of one of them 
(that corresponding to Ta{y, A^^)) start at the points y of 7 ,, and end at the points xi of the 
surface surrounded by 7 ^. The other bundle of parallel lines (corresponding to T^(x, A^^)) 
also begins at 7 *, but finishes at X 2 inside As before, the contributions to expression 
are expressed as areas of the infinitesimal parallelograms formed by vertical displacements 
along both the fiber and the curve where the corresponding fiber ends. 

Up to this point, there is nothing specially new in the discussion of the first-order invari¬ 
ant. The main difference with the previous case, is due to the presence of the loop coordinate 
with two indices which introduces an intrinsic order in curve i that affects which 

contributions must be summed up. In the integration along curve i the areas of the par¬ 
allelograms coming from curve k, which have to be multiplied by those coming from curve 
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j, contribute only if the former are “created” before the later when one travels along curve 
i starting at an arbitrary marked point {O, in FIG121). Henceforth, of the four products of 
areas |Hi||H 4 |, — |Hi||H 3 |, IH 2 IIH 3 I and — IH 2 IIH 4 I, that would appear if instead of 
we had times in the last one does not contribute due to the fact that A 2 is 
“drawn” after A 4 was. 

Then, the result of applying (|f)^ to the picture in FIG. El is minus the square of the area 
A‘f of any of the lobules of intersection. Glearly, this is an area-preserving diffeomorphism 
invariant, which could not be detected by the previous order one (which in fact vanishes). 

Summarizing, we have shown that classical 2d Yang-Mills theory coupled to non- 
dynamical Wong particles in Euclidean space can be used to obtain area-invariants through 
a perturbative scheme. This fact may be seen as an instance of a general result established 
in j^, which in turn aroused as a formalization of previous results jj, |^. 

The hrst two invariants provided by this method were explicitly computed, and a 
geometrically appealing interpretation was given. It is conjectured that the sequence of 
invariants here obtained could be seen as a kind of “projection” of the link-invariants of 
Milnor This and other pertinent questions will be addressed in future work. 


We thank Edmundo Gastillo for his assistance with the hgures. This work was partially 
supported by Fonacit grant G2001000712. 
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